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The stability of networks is greatly influenced by their degree distributions and in particular by 
their broadness. Networks with broader degree distributions are usually more robust to random 
failures but less robust to localized attacks. To better understand the effect of the broadness of the 
degree distribution we study two models in which the broadness is controlled and compare their 
robustness against localized attacks (LA) and random attacks (RA). We study analytically and by 
numerical simulations the cases where the degrees in the networks follow a bi-Poisson distribution 

P{k) = ae~^^ ^ ^ Gaussian distribution P(k) = A ■ exp(— ) 

with a normalization constant A where fc > 0. In the bi-Poisson distribution the broadness is 
controlled by the values of a, Ai, and A 2 , while in the Gaussian distribution it is controlled by the 
standard deviation, a. We find that only when q = 0 or a = 1, i.e., degrees obeying a pure Poisson 
distribution, are LA and RA the same. In all other cases networks are more vulnerable under LA 
than under RA. For a Gaussian distribution with an average degree fi fixed, we And that when is 
smaller than p the network is more vulnerable against random attack. However, when is larger 
than p the network becomes more vulnerable against localized attack. Similar qualitative results 
are also shown for interdependent networks. 


I. INTRODUCTION 

Complex networks are widely used as models to un¬ 
derstand such features of complex systems as structure, 
stability, and function [IMII. The robustness of networks 
suffering site or link attacks is a topic of great interest 
because it is an important issue affecting many real-world 
networks. Such approaches as site percolation on a net¬ 
work where nodes suffer either random attack (RA) 

14 or targeted attack (TA) based on node connectivity 

1^ have been developed to study these phenomena. 
Localized attack (LA) in which nodes surrounding a seed 
node are removed layer by layer has also been recently in¬ 
troduced [ 21 I, [ 2 ^ . In addition, interdependent networks 
are more vulnerable to RA and TA than isolated sin¬ 
gle networks [23l - l^ . LA on spatially embedded interde¬ 
pendent networks has been addressed, and a significant 
metastable regime where LA above a critical size propa- 
gcdes throughout the whole system has also been found 

Although prior research has developed tools for prob¬ 
ing network robustness against all these attack scenar¬ 
ios and has found that degree distribution broadness 
strongly influences network stability 0, there has been 
no systematic study of how degree distribution broad¬ 
ness affects robustness. Here we compare LA and RA on 
two networks models in which the broadness is controlled. 
One model is bi-Poisson with two groups having different 
average degrees. The difference between the two average 
degrees characterizes the broadness of the degree distri¬ 
bution of the network. Although research on this topic 
usually focuses on a network with a pure Poisson degree 
distribution, many real-world networks have two or more 
degree distributions HIMl. For example, a network of 
two groups of people, a high-degree group with many 


friends and a low-degree group with few friends, might 
reflect a bi-Poisson distribution. Note that bi-Piossonian 
networks are optimally robust against TA [s^. The sec¬ 
ond model in which the broadness can be controlled is a 
Gaussian degree distribution. Here the standard devia¬ 
tion a characterizes the broadness of the degree distribu¬ 
tion. This distribution is realistic, e.g., the distribution 
of WWW links resembles a Gaussian distribution . 

We here analyze the robustness against attack of net¬ 
works in which we can tune the broadness of the degree 
distributions, e.g., those with bi-Poisson and Gaussian 
degree distributions. We limit our approach to LA and 
RA and use the frameworks developed in Refs. 0 and 
[2l[, extending them to study (i) single networks with a 
bi-Poisson distribution, (ii) single networks with a Gaus¬ 
sian distribution, (iii) fully interdependent networks with 
the same bi-Poisson distribution in each network, and 
(iv) fully interdependent networks with the same Gaus¬ 
sian distribution in each network. By changing a of the 
bi-Poisson distribution 

P{k) = ae-^^ ^ + (1 - € P, 1], (1) 

with fixed Ai and A 2 , and of the Gaussian distribution, 

P(fc)=H.exp(-^^^^),fc>0, (2) 

with p fixed, we investigate how the distribution broad¬ 
ness influences the percolation properties. These include 
the size of the giant component Poo as a function p, the 
fraction of unremoved nodes and the critical threshold Pc 
at which the giant component Poo first collapses. In all 
cases we find that our extensive simulations and analyt¬ 
ical calculations are in agreement, and observe the qual¬ 
itative characteristics of robustness in both single and 
interdependent networks under both LA and RA. 
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II. RA AND LA ON A SINGLE NETWORK For a Gaussian distribution we have 


A. Theory 

Following Ref. [s^, we introduce the generating func¬ 
tion of the degree distribution P{k) of a certain network 
as 

Go(x) = ^P(fc)a;'=. (3) 

k 


P- = k\k - l)e(-G-/^)V2<.^) • 

The size of the resultant giant component is Q 

Pooip) = 1 - i?o(l) = p[l - Go(i?i(l))], (13) 

which can be numerically determined by solving -ffi(l) 
from its self-consistent equation 


Similarly, for the generating function of the underlying 
branching processes, we have 


Gi(x)=^ 

k 


P{k)k 

(k) 


G'oix) 

Goil)' 


(4) 


The size distribution of the clusters that can be reached 
from a randomly chosen link is generated in a self- 
consistent equation 


Hi{x) = xGi{Hi{x)). (5) 


Then the size distribution of the clusters that can be 
traversed by randomly following a starting vertex is gen¬ 
erated by 


Ho{x) = xGo {Hi{x)). (6) 

Next we distinguish between random attack and localized 
attack. 

(I) Random Attack: An initial attack with the ran¬ 
dom removal of a fraction 1 — p of nodes from the net¬ 
work changes the cluster size distribution of the remain¬ 
ing network and the generating functions of the surviving 
clusters’ size distribution become 0 

Hi{x) = l-p + pxGi{Hi{x)), (7) 

and analogously, 

Ho{x) = l-p + pxGo{Hi{x)). (8) 

Here Pc, the critical value at which the giant component 
collapses, is determined by 


Hi(l) = l-p + pGi(Hi(l)). (14) 


(II) Localized Attack: We next consider the local re¬ 
moval of a fraction 1 — p of nodes, starting with a ran¬ 
domly chosen seed node. Here we remove the seed node 
and its nearest neighbors, next-nearest neighbors, next- 
next-nearest neighbors, and continue until a fraction 1—p 
of nodes have been removed from the network. This 
pattern of attack reflects such real-world cases as earth¬ 
quakes or the use of weapons of mass destruction. As 
in Ref. [2l|, the localized attack occurs in two stages, (i) 
nodes belonging to the attacked area (the seed node and 
the layers surrounding it) are removed but the links con¬ 
necting them to the remaining nodes of the network are 
left in place, but then (ii) these links are also removed. 
Following the method introduced in Refs. [2l|,[U, we find 
the generating function of the degree distribution of the 
remaining network to be 


where / = Gg ^ (p). The generating function of the un¬ 
derlying branching process is thus 


Gpi{x) 



(16) 


The generating function of the cluster size distribution 
following a random starting node in the remaining net¬ 
work is 


Hpo{x)=xGpo{Hpi{x)), (17) 


and 


1 

1 g;(i) 

G;( 1 ) Gg(l)’ 


(9) 


( 10 ) 


which is equivalent to the expression pc = {k)/{k{k — 1)) 
given in Ref. Q- 

Thus for a bi-Poisson distribution, because Go(a;) = 

Q,gAi(x-i) _|_ — a)e^'^G-A ^ is 


aAi -l- (1 — q:)A2 

aX\ -I- (1 — q;)A2 


( 11 ) 


where Hpi(x), the generating function of the cluster size 
distribution given by randomly traversing a link, satisfies 
the self-consistent condition 


Hpi{x) = xGpi{Hpi{x)). (18) 

The network begins to generate a giant component when 
Gpi(l) = 1 [2l|, which yields Pc as the solution to 

Go(Go'(Pc)) = G;(1). (19) 

The size of the giant component Pooijp) as a fraction of 

the remaining network thus satisfies 

Poo(p)=p[l-Gpo(i7pi(l))], (20) 














FIG. 1. (Color online) Sizes of giant component, Poo(p), as 
a function of p for Ai = 4, A 2 = 12 and a — 0.7. Here 
solid lines are theoretical predictions, from Eq. (I13II for RA 
(red line) and Eq. (I20II for LA (green line), and symbols are 
simulation results with network size N = 10^, where averages 
are taken over 10 realizations, under LA (Q) and RA (□). 


which can be numerically determined by hrst solving 
Hpi{l) from Eq. (HH]), i.e., i?pi(l) = Gpi(i?pi(l)). 

In order to determine Pc explicitly, we first get fc from 
fc = Gq^{Pc), i-e., fc from Go(/c) = Pc-^ Then from 
Eq. (fTOl) fc must also satisfy Gq(/c) = Gg(l). In the 
general case, Pc and Poo must be obtained by solving 
numerically Eqs. (HU and (I20|) . In certain limiting cases, 
however, one can derive explicit analytical expressions 
for Pc that yield more physical insight. An example of a 
specihc case is given in the next subsection. 

1. Analytic solution of pc for bi-Poisson distribution with 
A 2 = 2Ai 

For a bi-Poisson distribution, using its generating func¬ 
tion and Go(/c) = Pc, fc and Pc satisfy the relation 

Goifc) = + (1 - = Pc- (21) 

Assuming A 2 = 2Ai, we denote — y |;]ga,t 

Eq. (1^ reduces to ay-\-{l — a)y^ = Pc, which, for a 1, 
is a quadratic equation of y and its positive solution is 

+ 4p.(l - (22) 

^ 2(1 - a) ^ ^ 

Plugging fc into Eq. ((11]) we get another quadratic equa¬ 
tion of y, 

cxX^^y -(- (1 — — cyXi -t- (1 — 

for which the physical solution of y is 

_ \/ a^Xf 4(1 — Q!)A|[a(Ai — A 2 ) -|- A 2 ] — aAf 
^ “ 2(l-a)Ai • 


FIG. 2. (Color online) Percolation thresholds pc of a single 
bi-Poisson network as a function of a under LA and RA with 
Ai = 4, A 2 = 12. Here solid lines are theoretical predictions, 
from Eq. HU for RA (blue line) and Eq. (fTOll for LA (green 
line) and symbols (□ for RA and Q for LA) with error bars 
are simulation results with network size of A = 10^ nodes, 
where averages and standard deviations are taken over 20 
realizations. 


Because fc = ln{y)/Xi-\-l, to obtainpc we need to equate 
Eqs. (1^ and (1^ . Thus we obtain 

Pc = g 4 ^ 2 ^_^^ [/^ + 6a\/o;^-k/3-6a^], (25) 

where f3 = ■ We use the relation of A 2 = 2Ai 

for simplification. Plugging a = 0 into Eq. (1251) . we get 
Pc = I/A 2 as found in Ref. [2lj. For a —>• 1, employing 
the L’Hopital rule we also get limQ_>i pc = 1/Ai, as found 
in the pure Poisson distribution described above. 

It is impossible to derive Pc explicitly for a Gaussian 
distribution. Even for a bi-Poisson distribution, other 
than special cases such as the one discussed above, de¬ 
riving Pc is also impossible because it requires solving 
first fc = Gq ^(pc), i.e., fc from Eq. (1^ . which could be 
viewed as ay^^ + (1 ~ = Pc, a polynomial equation 

of y = Because we also consider the cases of 

A 2 > Ai ^ 4 using the Abel-Ruffini theorem, there is no 
general algebraic solution to the above equation except 
in some special cases. Hence we use the Newton method 
to solve Pc and Poo numerically. 


B. Results 

To test the analytical predictions above we conduct nu¬ 
merical calculations of analytic expressions, and compare 
the results with the simulation results on single networks 
with degrees following both bi-Poisson distributions and 
Gaussian distributions under both LA and RA. All the 
simulation results are obtained for networks of A = 10^ 
nodes. 

















FIG. 3. (Color online) Comparison between numerical results 
(symbols) and the analytic results (solid lines) for bi-Poisson 
distribution with Ai = 4 and A 2 = 8. Note that they agree 
with each other well. Here, all the analytic results are ob¬ 
tained from Eq. (1111) for RA (red line) and Eq. (I25II for LA 
(black line) and the numerical results are attained by employ¬ 
ing Newton’s Method on Eqs. m and (O respectively. 



FIG. 4. (Color online) Sizes of giant component as a function 
of p of a single Gaussian network with p = 4 and = 2. 
Here solid lines are theoretical results, from Eq. (1131) for RA 
(blue line) and Eq. (I20II for LA (green line) and symbols are 
simulation results obtained from network size of A" = 10^ 
where averages are taken over 10 realizations under LA (Q) 
and RA (□). 


1. Single bi-Poisson networks 


Figure [T] shows the giant component Poo{p) as a func¬ 
tion of the occupation probability p under LA and RA. 
Note that Pc is larger for LA than for RA. The simu¬ 
lation results agree with the theoretical results obtained 
from Eqs. (USD and dSOD, and there is second-order perco¬ 
lation transition behavior in both attack scenarios. Note 
that when a = 0 or 1 , i.e., when node degrees follow 
a pure Poisson distribution as reported in Ref. [2l|, the 
networks have the same critical value of Pc under LA and 
RA and the same dependence of Poo{p) on p. However 
when a = 0.7, Pc{LA) > pc{RA), indicating that the net¬ 
work is more fragile under LA than under RA, and that 
the giant components exhibit different behavior. 

Figure [2] shows how the broadness of the distribution, 
tuned by changing a with fixed Ai and A 2 , influences the 
robustness of the network under LA and RA. The solid 
lines are the numerical results obtained from the Newton 
method and the symbols with error bars are the simula¬ 
tion results. Note that only when a = 0 and a = 1 does 
Pc{LA) = pc{RA). In all other cases pc{LA) > pc{RA), 
indicating that the network is always more vulnerable 
under LA than under RA if the degree distribution is 
bi-Poissonian. Note also that pc{LA) peaks at a = 0.79. 

For the special case of A 2 = 2Ai, we compare the 
analytical values of Pc from Eqs. (ED and (EID using 
Ai = 4 and A 2 = 8 with results obtained from the New¬ 
ton method (see Fig. |3D. For this combination of average 
degrees, pc{LA) peaks at a = 0.91. Note that the results 
agree, indicating that the Newton method produces satis¬ 
factory results and therefore, in the general case in which 
A 2 7 ^ 2Ai and in the cases of Gaussian distribution, it can 
be used to get pc{LA). 


2. Single Gaussian Networks 


Figure ID shows the giant component Poo{p) as a func¬ 
tion of the occupation probability p under LA and RA 
respectively for a single network with a Gaussian de¬ 
gree distribution. Note that the simulation results and 
the theoretical results obtained from Eqs. ED and (EOD 
agree, and that second-order phase transition behavior is 
present in both attack scenarios. Note also that p = 4 
and cr^ = 2 , and thus Pc{LA) < pc{RA), which indicates 
that the network is more robust under LA than under 
RA for this particular distribution. 

We fix p and find that when the Gaussian distribution 
gets broader, i.e., when a increases, Pc{RA) decreases, 
but that Pc{LA) increases with a (see Fig. [SD. Note that 
when tr^ < p, pc{LA) < pc{RA), and that the opposite 
is true when cr^ > p. Note also that when p there 
is a crossing point with pc{RA) « pc{LA), which is anal¬ 
ogous to a Poisson ER network with the same mean and 
variance and the robustness of the network under both 
LA and RA is the same, as reported in Ref. plj . 

Figure [5] shows a plot of as a function of p when this 
intersection point occurs, i.e., when pc{LA) = pc{RA). 
Note that except for some minor deviations at small p 
values, because k > 0 the Gaussian distribution is de¬ 
formed, the region above the extrapolation curve corre¬ 
sponds to Pc{LA) > pc{RA), and the region below to 
Pc{LA) < Pc{RA). 










FIG. 5. (Color online) Percolation thresholds pc as a function 
of of networks with Gaussian degree distribution under LA 
and RA with /r = 4. Here solid lines are theoretical predic¬ 
tions, from Eq. (1121) for RA (red line) and Eq. (1191) for LA 
(black line) and symbols (□ for RA and Q for LA) with er¬ 
ror bars are simulation results with network size of = 10^ 
nodes, where averages and standard deviations are taken over 
20 realizations. It is shown here that as increases pc(LA) 



FIG. 6. (Color online) as a function of p, at the intersec¬ 
tion point where pc{LA) = pc{RA) for single networks where 
degrees follow a Gaussian distribution. 


III. RA AND LA ON FULLY 
INTERDEPENDENT NETWORKS 

A. Theory 

We apply the formalism of RA on fully interdependent 
networks introduced in Ref. [2^ . Specifically, we consider 
two networks A and B with the same number of nodes N. 
Within each network the nodes are randomly connected 
with degree distributions PA{k) and Psik) respectively. 



FIG. 7. (Color online) Sizes of the mutually connected giant 
component of the fully interdependent bi-Poisson networks 
as a function of p for Ai = 4, A 2 = 12 and a = 0.5. Here 
solid lines are theoretical predictions, from Eq. (1261) for RA 
(blue line) and similarly for LA (green line), and symbols are 
simulation results with network size N = 10"^, where averages 
are taken over 10 realizations, under LA (Q) and RA (□). 

Every node in network A depends on a random node in 
network B, and vice versa. We also assume that if a 
node i in network A depends on a node j in network B 
and node j depends on node I in network A, then I = i, 
which rules out the feedback condition [^. This full 
interdependency means that every node i in network A 
has a dependent node j in network B, and if node i fails 
node j will also fail, and vice versa. 

(I) Random Attack: We begin by randomly removing 
a fraction 1 — p of nodes and their links in network A. 
All the nodes in network B that are dependent on the re¬ 
moved nodes in network A are also removed along with 
their connectivity links. As nodes and links are sequen¬ 
tially removed, each network begin to break down into 
connected components. Due to interdependency, the re¬ 
moval process iterates back and forth between the two 
networks until they fragment completely or produce a 
mutually connected giant component with no further dis¬ 
integration. As in Ref. [ 2 ^ we introduce the function 
5a(p) = 1 — Gao[1 — p(l — /a)], which is the fraction of 
nodes that belong to the giant component of network A, 
where /a is a function of p that satisfies the transcenden¬ 
tal equation /a = Gai[1 — p(l — /a)]- Similar equations 
exist for network B. When the system of interdependent 
networks stops disintegrating, the fraction of nodes in the 
mutual giant component is Poo, satisfying 

Poo = xqb^x) = ygA{y), (26) 

where x and y satisfy 

x = pgA{y),y =pgB{x). (27) 

Excluding the trivial solution a; = 0, p = 0 to the equa¬ 
tion set above, we combine them into a single equation 











FIG. 8. (Color online) Percolation thresholds Pc of the fully 
interdependent bi-Poisson networks with Ai = 4, A 2 = 12 
as a function of a under LA and RA. Here solid lines are 
theoretical predictions, from Eq. (1291) for RA (bine line) and 
similarly for LA (green line) and symbols (□ for RA and Q for 
LA) with error bars are simulation results with network size 
of A = nodes, where averages and standard deviations 
are taken from 20 realizations. When a is not 1 or 0, pc{LA) 
is always larger than pc{RA). 


by substitution and obtain, 


X = gA[9B{x)p]p. (28) 


A nontrivial solution emerges in the critical case (x = 
Xc,p = Pc) by equating the derivatives of both sides of 
Eq. (|^ with respect to x 


1 =p- 


dgAlpgsix)] dgsix) 


dx 


dx 


X=Xa,p=Pc 


(29) 


which, together with Eq. (123, gives the solution for pc 
and the critical size of the mutually connected compo¬ 
nent, Poo{Pc) = Xcgsixc). 

(II) Localized Attack: When LA is performed on the 
one-to-one fully interdependent networks A and B de¬ 
scribed above, we can hnd an equivalent random net¬ 
work E with generating function Ge{){x) such that after 
a random attack in which 1 —p nodes in network E are re¬ 
moved, the generating function of the degree distribution 
of the remaining network is the same as Gpo ( 2 ;) (with the 
substitution of Go(a;) by Gaq{x)). Then the LA problem 
on networks A and B can be mapped to a RA problem on 
networks E and B. By using G£;o(l — P + px) = Gpo{x) 
and from Eq. ca we have 


G_Eo(a;) 


Gao(/) 


G^o[/- 


GaoU) 


G'^^{l)GA^{f) 


(x-1)]. (30) 


Thus by mapping the LA problem on interdependent 
networks A and R to a RA problem on a transformed pair 
of interdependent networks E and B, we can apply the 
mechanism of RA on interdependent networks to solve Pc 
and Poo(p) under LA. 


FIG. 9. (Color online) Sizes of mutual giant component of the 
fully interdependent Gaussian networks as a function of p with 
p = 4 and cr^ = 2. Here solid lines are theoretical predictions, 
from Eq. (12611 for RA (red line) and similarly for LA (black 
line), and symbols are simulation results with network size 
N = 10“^, where averages are taken over 10 realizations, under 
LA (O) and RA (□). 


Note that for pure Poisson distributions, / = 

G^lip) = and that by substituting / into 

Eq. (1501) we get Geo{x) = Gao{x). Thus we find that 
pure Poisson distributions have exactly the same perco¬ 
lation properties for fully interdependent networks under 
LA as those under RA, as found in Ref. [ 2 ^. Because 
the extreme complexity of the above equations makes it 
difficult to obtain explicit expressions for pc and Poc{p) 
except when degree distributions are simple, we resort to 
numerical calculations in general. 


B. Results 

1. Fully interdependent networks with bi-Poisson degree 
distribution 

We start with two fully interdependent networks in 
which the degrees both follow the same bi-Poisson distri¬ 
bution and carry out a RA on one of the networks, initi¬ 
ating a cascading failure process that will continue until 
equilibrium is reached. We then do the same procedure 
with the same set-up but this time using a LA to initiate 
the cascading failure process. Figure [7] shows the size of 
the giant component Poo{p) of the system as a function 
of the occupation probability p under LA and under RA. 
Note that in both RA and LA scenarios the simulation re¬ 
sults and the theoretical results obtained from Eq. (1551) 
agree, indicating that our strategy of finding an equiv¬ 
alent network under LA works. The first-order phase 
transition that occurs in both attack scenarios indicates 
that the interdependency of the system makes it much 
more vulnerable to attack than single networks. When 
a = 0.5 the system is more fragile under LA than under 
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FIG. 10. (Color online) Percolation thresholds pc as a func¬ 
tion of of the fully interdependent Gaussian networks un¬ 
der LA and RA with /r = 4. Here solid lines are theoretical 
predictions, from Eq. (1291) for RA (red line) and similarly 
for LA (black line) and symbols (□ for RA and O for LA) 
with error bars are simulation results with network size of 
N = 10^ nodes, where averages and standard deviations are 
taken from 20 realizations. It is seen here that as increases 
Pc{LA) increases and Pc{RA) has a tendency to decrease. As 
cr^ approaches the value of p, pc{LA) « pc{RA), which is 
manifested by the intersection point shown here. 


RA with pc{LA) > pc{RA), and the giant components 
exhibit different behaviors. 


Figure [5] shows how the broadness of the distribu¬ 
tion, tuned by changing a with fixed Ai and A 2 , influ¬ 
ences the robustness of the network under both LA and 
RA. Solid lines are numerical results using the Newton 
method on Eq. (1^51) and symbols with error bars are 
simulation results. Note that only when a = 0 and 
a = 1 is P{k) reduced to a pure Poisson, and we have 
Pc{LA) = pc{RA) = 2AbbAl (k), as in Ref. [ 2 ^. When 
a deviates from 0 or 1, i.e., when P{k) deviates from 
a pure Poisson distribution and takes the form of a bi- 
Poisson distribution, pc{LA) > pc{RA), indicating that 
the system is more vulnerable under LA than under RA. 


2. Fully interdependent networks with Gaussian degree 
distribution 

We construct two fully interdependent networks in 
which the degrees in each network follow the same Gaus¬ 
sian distribution and carry out a RA on one of the net¬ 
works to initiate a cascading failure process that will con¬ 
tinue until it reaches a steady state. We repeat the ac¬ 
tion, but this time using a LA. Figure IHl shows the sizes 
of the giant component Poo (p) as a function of the occu¬ 
pation probability p under both LA and RA. Note that 
simulation results and the theoretical results obtained 
from Eq. (E51) agree. When /i = 4 and cr^ = 2 the 
system is more fragile under LA than under RA with 
Pc{LA) < pc{RA), and the giant components exhibit dif¬ 
ferent behaviors. 

If we fix p, when the Gaussian distribution gets 
broader, i.e., when a increases, analogous to what we 
find in a single Gaussian network, the critical Pc be¬ 
havior of the system differs under LA from that un¬ 
der RA. Figure [TO] shows the effect of cr on pc in the 
fully interdependent Gaussian networks. When cr^ < p, 
Pc{LA) < pc{RA), and the opposite occurs when > p. 
The intersection point in Fig. [TUlis located near p, 
similar to that in Poisson distribution networks. Thus 
the system behaves the same under LA as under RA, con¬ 
firming the results presented in the previous subsection. 
Note that our results show that in both attack scenarios, 
the interdependency of the system makes it much more 
vulnerable to RA and LA compared to single networks 
(compare Fig. [TUlto Fig. [5]). 

IV. CONCLUSIONS 

In summary, we show that a LA on interdependent net¬ 
works can be mapped to a RA problem by transforming 
the network under initial attack. We also show how the 
broadness of the degree distribution affects the robust¬ 
ness of networks against RA and LA respectively. We 
show that, in general, as the degree distribution broad¬ 
ens the network becomes more vulnerable to LA than 
RA. This finding holds for both single networks and in¬ 
terdependent networks. 
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